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Abstract 

We consider nematic liquid crystals in a bounded, convex polyhe- 
dron described by a director field n(r) subject to tangent boundary 
conditions. We derive lower bounds for the one-constant elastic en- 
ergy in terms of topological invariants. For a right rectangular prism 
and a large class of topologies, we derive upper bounds by introducing 
test configurations constructed from local conformal solutions of the 
Euler-Lagrange equation. The ratio of the upper and lower bounds 
depends only on the aspect ratios of the prism. As the aspect ra- 
tio is varied, the minimum-energy conformal state undergoes a sharp 
transition from being smooth to having singularities on the edges. 
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The continuum theory of liquid crystals [l a is a prototypical nonlinear field 
theory in which topological considerations play a fundamental role El H] , 
both in equilibrium (eg [SHE]) and dynamical phenomena (eg [Jj). Nematic 
liquid crystals are represented by a director field n(r), which describes the 
mean local orientation of the constituent rod-like molecules. In confined 
geometries, boundary conditions, which depend on the properties of the sub- 
strate, can play a significant role (see, eg ||). 

In this Letter, we consider director fields in a bounded, three-dimensional 
convex polyhedron P. While natural theoretically, the problem also has 
a technological motivation; polyhedral geometries have been proposed as 
a mechanism for engendering bistability in liquid crystal display cells [§], 
[TU] . Polyhedral cells can support two (or more) energetically stable director 
configurations with contrasting optical properties. Power is required only to 
switch pixel states but not to maintain them. Bistable cells offer the promise 
of displays with higher resolution and requiring less power than is available 
with current technologies based on monostable cells. 

We suppose P has strong azimuthal anchoring, so that n satisfies tangent 
boundary conditions - on the faces of P, n must be tangent to the faces. 
Tangent boundary conditions imply that, on the edges of P, n is parallel to 
the edges, and, therefore, is necessarily discontinuous at the vertices. We 
restrict our study to configurations which are continuous everywhere else, ie 
as continuous as possible. As P is simply connected, we can then regard 
n(r) as a unit-vector field, rather than a director field. We obtain lower 
bounds for the elastic energy of n in terms of its topological invariants, and, 
for a rectangular prism, upper bounds, which differ from the lower bounds 
only by a geometry-dependent factor. The upper bounds are obtained from 
local conformal solutions of the Euler-Lagrange equation, whose energetics 
indicate the onset of edge singularities as the prism becomes cubic. 

It turns out that tangent boundary conditions produce a large family 
of topologically distinct configurations. A complete classification is given in 
Ref. whose results we briefly summarise. Tangent unit- vector fields on P 
can be classified up to homotopy (ie, continuous deformations) by a family 
of invariants: the edge orientations, kink numbers and trapped areas. The 
edge orientations are just the values of n on the edges of P, and therefore are 
essentially signs. The kink numbers determine the number of times n winds 
along an (outward-oriented) path on a face of P which joins a pair of adjacent 
edges. On such a path, the initial and final values of n are determined by 
the edge orientations; in between, n describes a curve on the circle of unit 
vectors tangent to the face. The shortest curve between the endpoints is 
assigned kink number 0. In general, the kink numbers are integers. The 
trapped areas, fl a , are defined as follows. Let C a be a surface inside P which 
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separates the vertex a from the other vertices. Then tt a is the area of the 
region n(C a ) on the unit two-sphere. This may be written as j ca D ■ dS, 
where dS is the outward-oriented area element on C a , and 

D 3 = ^jki(d k n x din) • n. (1) 

D may be regarded as the vector field dual to the pull-back, n*uj, of the 
outward-oriented area form u on the sphere (in polar coordinates, ui = 
sin 6 d9 A d<f>). fl a is independent of the choice of C a , and may be thought of 
as the degree of a fractional point defect at a. fl a need not be integral, but 
the allowed values of Q a for given edge orientations and kink numbers differ 
by integer multiples of 47r (complete coverings of the sphere). These integers 
are the wrapping numbers. The invariants satisfy certain sum rules, which 
follow from the fact that n is continuous away from the vertices. The sum of 
the kink numbers on each face is determined by the edge orientations, while 
the trapped areas sum to zero. All values of the invariants consistent with 
the sum rules can be realised, and two configurations are homotopic if and 
only if their invariants are the same. 

In the continuum theory, the energy of n is given by the elastic, or Frank- 
Oseen, energy [I], 

E = J [i^i(divn) 2 + K 2 {n ■ curl n) 2 + K 3 (n x curl n) 2 

+ K 4 div((n- V)n- (divn)n)] dV. (2) 

Tangent boundary conditions imply that the contribution from the .fQ-term, 
which is a pure divergence, vanishes. In the so-called one-constant approxi- 
mation, the remaining elastic constants K\, K 2 and K 3 are taken to be the 
same. In this case, © simplifies to 




Stable (or quasi-stable) configurations are minima (or local minima) of the 
energy. In the one-constant approximation, local minimisers satisfy the 
Euler-Lagrange equation and boundary conditions 

An = (n • An)n, (4) 
(c • n) dP = (c x V c n) 9P = 0, 

where c is the outward unit-normal on the boundary dP of P, and V c denotes 
the derivative along c (the normal derivative here); the condition on V c n 
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ensures there is no boundary contribution to the first-order variation of the 
elastic energy. (Minimisers of the energy (jSJ) are called harmonic maps in the 
mathematics literature.) 

Our first result is a lower bound for the energy (J3J) in each homotopy class. 
The derivation extends an argument from a seminal paper of Brezis, Coron 
and Lieb [3] to the polyhedral boundary-value problem we are considering 
here. Below we restrict to the case where Vn is smooth away from the ver- 
tices. In fact, the bound is good provided only that Vn is square-integrable 

m 

Since n is normalised, the partial derivatives 9,-n, j = 1, 2, 3, are orthog- 
onal to n, and therefore are linearly dependent. Thus, at each point r, there 
is at least one direction, w say, in which the derivative of n vanishes. Let us 
introduce unit vectors u and v which together with w define an orthonormal 
frame. Then 

(Vn) 2 = (V u n) 2 + (V,n) 2 > 2|V u n x V„n|, (5) 

where equality holds if and only if V M n and V^n are orthogonal and of equal 
length. From the right-hand side of (jnj) is just 2|(u x v) • D|, which in 
turn is just equal to 2|D|, since, by inspection, D is orthogonal to u and v. 
Then 

(Vn) 2 > 2|D|. (6) 

Let £ (r) be a continuous scalar function with piecewise continuous gradi- 
ent such that |V£| = 1. Then |D| > D • V£. Let P be a domain obtained by 
excising from P vanishingly small neighborhoods of each vertex. From (JHJ) 
and ©, it follows that 

E > 2K j * D- Vf dV. (7) 

But D has zero divergence, as follows from the fact that area-form u is closed 
(as well as from direct calculation). Integrating by parts, we get 

E > 2K I £D • dS. (8) 

JdP 

The boundary of P, denoted dP, consists of (parts of) the faces of P as well 
as the boundaries of the excised regions. But D vanishes on the faces of P, 
as tangent boundary conditions imply that the tangential derivatives of n 
are parallel to each other there. On the boundaries of the excised regions, 
£ can be set to its values £ a at the vertices. Thus, the integral J d p^D ■ dS 
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reduces to a weighted sum of fluxes of D through the excised boundaries. 
These fluxes are just the trapped areas. We obtain the lower bound 



E > 2K^2CV a . 



(9) 



a 



To proceed, note the values £ a are constrained; since |V£| = 1, we must 
have |£ a — £ b \ < L ab , where L ab is the length of the edge between vertices 
a and b. Conversely, given a set of values £ a satisfying these inequalities, 
we can construct a function £(r) which interpolates between these values, 
and for which V£ is piecewise continuous with unit norm (for example, let 
£(r) = max a (£ a — |r — a|)). Therefore, we can express the lower bound 
in terms of the trapped and the lengths of the edges of P, 



(jlOjl specifies a linear optimisation problem, which can be solved by standard 
methods. 

We note that because Q a = 0, the optimal solution £ a is in general 
only determined up to an overall additive constant. If the one-constant ap- 
proximation is dropped, K in (jID}) can be replaced by the smallest of the 
three elastic constants Ki, K2 and K%. 

Next, we consider the case where P is a right rectangular prism, with sides 
of length L x > L y > L z . Moreover, within P, we consider configurations (and 
their homotopy classes) which are reflection- symmetric about the mid-planes 
of P. That is, 

n(L x - x, y, z) = n(x, L y -y,z) = n(x, y, L z - z) = n(x, y, z) (11) 

(we plan to treat more general geometries and homotopy classes in the fu- 
ture). Clearly, a reflection-symmetric configuration is determined by its val- 
ues in the octant of the prism, denoted R = {0 < Tj < \Lj}, as are its 
invariants. In particular, the trapped areas at two vertices are either the 
same or differ by a sign according to whether the vertices are related by an 
even or odd number of reflections. Let Q° denote the trapped area at the 
origin. Taking £ a in ([10)1 (optimally) to be L z or according to whether Q a 
is \Q°\ or — we obtain the explicit lower bound E > E_, where 



In the remainder of this Letter, we introduce a family of reflection- 
symmetric configurations which are local solutions of the Euler-Lagrange 




(10) 



= 8KL z \a 



(12) 
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equation Q. From these we infer upper bounds for the energy as well as the 
onset of edge singularities as the geometry is varied. 

Candidates for low-energy configurations are those for which (JHJ) becomes 
an equality. This is so if, on R, n is radially constant (ie, n(Ar) = n(r)) and 
conformal. Here, conformal means that the map t > V*n(r) from vectors t 
orthogonal to f to vectors V^n(r) orthogonal to n(r) preserves orientation, 
angles and ratios of lengths (or else vanishes). If n is radially constant then 
D is radial (cf (fl|). while if n is also conformal, D = |(Vn) 2 f. It is straight- 
forward to verify that radially constant, conformal configurations satisfy (J1J 
in R (but fail to on the mid-planes of P, where the normal derivatives of n 
are discontinuous). 

Suppose n is a reflection-symmetric configuration which is radially con- 
stant and conformal in R (below we will show how to construct such config- 
urations). Then, from (jH)), and using reflection symmetry, we get 

E = 8K I (Vn)W = 16K [ r • BdV = 16K [ rD • dS, (13) 

JR JR JdR 

where we have used V • D = 0. We obtain an upper bound E + for E by 
replacing r by its maximum value on dR, namely \{L 2 X + L 2 + Li) 1 / 2 . The 
integral which remains is just the flux of D through R. As tangent boundary 
conditions imply there is no contribution from the exterior faces of R (ie, 
x = 0, y = or z = 0), this just gives the trapped area. Thus 

E + = 8K(Ll + Ll + Ll) 1 ' 2 \rt>\. (14) 

Like E_, E + is proportional to and E + /E_ = (a 2 xz + a 2 z + l) 1 / 2 , where 
ctij = Li I Lj denote the aspect ratios of the prism. For the cube, this ratio is 

Vs. 

It remains to construct radially constant, conformal configurations n(r) 
in R satisfying tangent boundary conditions. For this, it is convenient to 
introduce the stereographic projection (e x , e y , e z ) \— > (e x + ie y )/(l + e z ) from 
the sphere to the extended complex plane C e . Projecting both r and n, we 
obtain a complex- valued function f(w) of complex argument w given by 

(15) 

The domain of / is the stereographic image of R, which is the quarter-disc Q 
given by |iw| < 1 and Rew, Imw > 0. n conformal is equivalent to / locally 
analytic. 

The form of / is determined by tangent boundary conditions. Under 
stereographic projection, the xz-iace is mapped to the real w-axis, so we 
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require that i) f(w) is real if w is real. Similarly, the yz- and xy-faces are 
mapped to the imaginary axis and the unit circle respectively, so we require 
that ii) f(w) is imaginary if w is imaginary, and iii) \ f(w) \ = 1 if |iu| = 1. To 
proceed, we assume that / has a meromorphic extension from Q to C e . Then 
the conditions i)— iii) may be continued to the following functional equations: 

f(w) = f{w), (16a) 
f(-w) = -f(w), (16b) 
f(w)f(l/w) = 1. (16c) 

Together, ()16a|) and ()16c|) imply that f(w)f(l/w) = 1, which in turn implies 
that if w is a zero of /, then 1/w is a pole. / meromorphic implies the 
number of zeros and poles in the unit disk must be finite (otherwise there 
would be an accumulation of one or the other). It follows that / has only a 
finite number of zeros and poles in C 6 , so that / is rational. 

(j!6b|) and ()16cj) imply that, if w is a zero of /, then so are w and — w, 
and similarly if w is a pole. One can then show that every reciprocal pair 
{w , 1/w) of zeros and poles has a unique representative in Q. For convenience, 
we divide these representatives into three groups: the strictly real r±, . . . , r a , 
with < Tj < 1 (as (|16cj) rules out zeros and poles of modulus one); the 
strictly imaginary is\, . . . ,isb, with < sj~ < 1; and the strictly complex 
ti,...,t c , with < \ti\ < 1. In addition, w = is a zero or pole of odd order, 
as is implied by ()16aj) and (j!6bjl . Then / is given by 

j— I \ J / f, = l \ K / 

LL{ {t 2 w 2_ 1)( fi w 2_ 1 

Here, pj is 1 or —1 according to whether rj is a zero or a pole (note that 
the Tj's needn't be distinct), and similarly for and 77. n, an odd integer, 
determines the order of the zero or pole at the origin. The coefficient e, given 
by (— l) a /(l), is equal to ±1 (from ()16aj) and (j!6c|l ). It is straightforward to 
verify that ()17|l satisfies tangent boundary conditions, and therefore charac- 
terises, via ()15|). the radially constant, conformal, tangent unit- vector fields 
n in R for which / has a meromorphic extension to C e . We remark that 
rational functions have been used extensively as the basis for an ansatz for 
skyrmions, minimum-energy configurations of fixed baryon number (ie, de- 
gree) of localised, nonlinear SU (2)- valued fields in M 3 . The ansatz yields very 
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good descriptions of the (numerically determined) energies and topologies of 
the true minimisers (see, eg, [T3], [H], [To]). 

The values of the topological invariants can be computed in terms of the 
parameters of /. From straightforward calculation gives D = Ar, where 
A, regarded function of w, is given by 

^ = 4in 2 /(i + i/i 2 ) 2 . as) 

It follows that Q° = JqAc1 2 w. The conditions (|16jl above imply that A is 
invariant under w i— > — w, w t— ► ty, and w i— > 1/iy. Therefore, 8f2° is given 
by the integral of ^4 over C e . But this quantity is just {An times) the degree 
of /, regarded (inverse stereographically) as a map from S 2 into itself. The 
degree of a meromorphic function is just the number of its zeros (or any other 
value) counted with multiplicity. Therefore, 

Q° = \{\n\ +2(a + b) +4c)tt. (19) 

The edge orientations are easily determined from the values of / at 1, i and 
(the images of the x, y and ,2-edges respectively). Let e x = ±1 according 
to whether n = ±x along the x-edge, and similarly for e y and e z . Then 

e x = e(-l) a , e y = e (-l) 6 (-l)('- 1 )/ 2 , e z = sgnn. (20) 

The kink number k z along the z-face can be computed from the change in 
phase of / along the quarter circle of |iu| = 1 between 1 and i. The kink 
numbers k x and k y along the x and y -faces can be computed by counting 
zeros of / with an appropriate sign along the respective intervals [0, 1] and 
[0, i\. The result is 

K = -|(-i)%(ELi(-i) fc ^ + |(i-(-i) 6 )^), 

ky = -|(-l)°e £C (E-= 1 (-l)^ + 1(1 "(-!)>,), (21) 
k z = \ {e x e y — n) — \Yl,j=iPj ~ \Hk=i a k — J2i=i T i- 

As remarked above, the kink numbers and edge orientations determine the 
trapped areas up to an integer multiple of Att. Moreover, trapped areas 
with any such multiple can be realised. For conformal configurations (fTTj) . 
it can be shown, using (|19j) - (J21j) . that all values of the trapped area fl° 
greater than f2 min = 27r(|fc. r | + \ k y \ + \k z \ + |) can be realised. Anticonformal 
configurations (obtained by replacing w with w) allow for all values of Q° less 
than — f2 min . Further details along with a discussion of the remaining values 
of fl° will be given elsewhere [To] . 

It is possible to get better upper bounds by evaluating the integrals in (fTH|) 
exactly. We do this explicitly for the configuration f(w) = w. We call this 



8 



configuration unwrapped, because it is topologically the simplest possible; 
its kink numbers are zero and its trapped area has the minimum allowed 
value. There are in fact eight distinct unwrapped configurations generated 
by the transformations w i— > —w, 1/w and w. All have the same energy. For 
f(w) = w, it is readily computed that D = r/r 3 . From (j!3)l . the contribution 
to the energy from the interior face {x = \L X \ of R is given by 



8KL X 



r 2 Lz 



dydz 



i i 



= 2a yx a zx KL x / / — — ■ — - — du dv. (22) 

Jo Jo ± + a 2 yx u + a 2 zx v 

The last integral may be identified with an integral representation of the 
Appell hypergeometric function F 2 (a, (3, /3', 7, 7'; s, t) [T^j, with parameters 
a = 1, (3 = ft' = ~, 7 = 7' = I and arguments s = —a 2 x and t = —a 2 zx . The 
integrals over the two other interior faces of R are evaluated similarly, while 
the integrals over the external faces vanish by tangent boundary conditions. 
Thus, the unwrapped energy is 

-^0 = 8 2^ a ji a kiKL i F2{l- l \, I , §, §, — OfcJ, (23) 

i 

where % runs over x, y and z, and k) is a cyclic ordering of (x, y, z). For 
a cube ( L = K — 1), this gives an upper bound of 15.3, about 20% more 
than E- = An. 

For conformal configurations other than the unwrapped ones, the energy 
16K f dR r~D ■ dS depends on the positions rj, Sk, t t of the zeros and poles, 
and therefore can be minimised with respect to these parameters. Since 
| aj jD • dS is the trapped area Q°, and therefore is parameter- independent, 
it is evident that the minimum energy is achieved by making |D| small at 
points of dR where r is large, and large at points where r is small. The 
local minima of r are just the corners of R on the x-, y- and z-edges, with 
projections w — 1, i and 0. For a cube, the corners are all equally close 
to the origin, and the minimum energy is approached in the singular limit 
in which the zeros and poles of / are made to coalesce, pairwise, at 1, % or 
(so that, from ()18|1 . |D| = A is made to diverge there), leaving a single 
zero or pole at w = 0. In this limit, all the topologically nontrivial behaviour 
(kinking and wrapping) concentrates at the edges, while away from the edges 
the configuration becomes unwrapped. This is reminiscent of the dipole con- 
figurations of Ref . [H] . However, for rectangular prisms, the minimum energy 
within a conformal family may be realised for a nonsingular configuration. 
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Figure 1: Scaled energy e = E/V 
Solid curve: L~ - 



s 2 )/(s 2 w 2 + 1 



of the conformal configuration w(w 2 + 
20, L y = 10, L z = 1. Dashed curve: 



This is illustrated in Fig. ^ which shows the energy (scaled by the cube 
root of the volume) of the configuration f(w) = w{w 2 + s 2 )/(s 2 w 2 + 1). For 
a cube (dashed curve), the energy approaches a minimum as s approaches 
1, corresponding to a configuration which is singular along the y-edge. For 
L x = 20, L y = 10, L z — 1 (solid curve), the energy has a minimum for s 
between and 1, corresponding to a smooth configuration. (The nontrivial 
kink number \k z \ = 1 rules out conformal deformations of / to a configura- 
tion which is singular along the shortest z-edge.) Numerical solutions of the 
Euler-Lagrange equations exhibit the same transition ^Sj. One would 
expect to observe the smooth configurations but not necessarily the singular 
ones; in the vicinity of such an edge singularity, the liquid crystal could melt, 
losing orientational order, and then relax to an unwrapped state. 

The preceding analysis indicates that the stability of topologically non- 
trivial tangent director configurations in rectangular prisms depends on the 
(purely geometrical) aspect ratios, a phenomenon which will be the subject 
of further investigation. 
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